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ABSTRACT The radiation we observe from blazars is most likely the product of the transforma- 
tion of bulk kinetic energy into random energy. This process must have a relatively small efficiency 
(e.g. 10%) if jets are to power the extended radio-structures. Recent results suggest that the av- 
erage power reaching the extended radio regions and lobes is of the same order of that produced 
by accretion and illuminating the emission line clouds. Most of the radiative power is produced in 
a well localized region of the jet, and, at least during flares, is mainly emitted in the 7-ray band. 
A possible scenario qualitatively accounting for these facts is the internal shock model, in which 
the central engine produces a relativistic plasma flow in an intermittent way. 
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1. INTRODUCTION 

We believe that the continuum radiation we see from blazars comes from the trans- 
formation of bulk kinetic energy, and possibly Poynting flux, into random energy 
of particles, which quickly produce beamed emission through the synchrotron and 
the inverse Compton process. This is analogous to what we believe is happening in 
gamma-ray bursts, although the bulk Lorentz factor of their flow is initially larger. 

Evidences for bulk motion in blazars with Lorentz factors between 5 and 20 
have been accumulated along the years, especially through the monitoring of super- 
luminally moving blobs on the VLBI scale (Vermeulen & Cohen 1994), and, more 
recently, through the detection of very large variable powers emitted above 100 MeV 
(see the third EGRET catalogue, Hartman et al., 1999), which require beaming for 
the source to be transparent to photon-photon absorption (e.g. Dondi & Ghisellini, 
1995). 

The explanation of intraday variations of the radio flux, leading to brightness 
temperatures in excess of T-q = 10 18 K (much exceeding the Compton limit) are 
instead still controversial (Wagner & Witzel 1995). Interstellar scintillation is surely 
involved, but it can work only if the angular diameter of the variable sources is so 
small to nevertheless lead to Tb = 10 15 K, which requires either a coherent process 
to be at work (e.g. Benford & Lesch 1998) or a Doppler factor of the order of a 
thousand. 
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Another controversial issue is the matter content of jets. We still do not know if 
they are dominated by electron-positron pairs or by normal electron-proton plasma 
(see the reviews by Celotti, 1997, 1998). 

Part of our ignorance comes from the difficulty of estimating intrinsic quantities, 
such as the magnetic field and the particle densities, using the observed flux, which 
is strongly modified by the effects of relativistic aberration, time contraction and 
blucshift, all dependent on the unknown plasma bulk velocity and viewing angle. 
Furthermore it is now clear (especially thanks to multiwavelength campaigns) that 
the blazar phenomenon is complex. 

On the optimistic side, we have for the first time a complete information of 
the blazar energy output, after the discovery of their 7-ray emission, and some 
hints on the acceleration process, through the behaviour of flux variability detected 
simultaneously in different bands (see the review by Ulrich, Maraschi & Urry 1996). 
Also, blazar research can now take advantage of the explosion of studies regarding 
gamma-ray bursts, which face the same problem of how to transform ordered to 
random energy to produce beamed radiation (for reviews: Piran 1999; Meszaros 
1999). 

2. ACCRETION = ROTATION? 

Despite the prediction that jets carry plasma in relativistic motion dates back to 
1966 (Rees, 1966), and intense studies over the last 20 years (Begelman, Blandford 
& Rees, 1984), quantitative estimates of the amount of power transported in jets 
have been done only relatively recently, following new observational results. 

One important point is that the extended (or lobe) radio emission of radiogalax- 
ies and quasars traces the energy content of the emitting region. Through minimum 
energy arguments and estimates of the lobe lifetime by spectral aging of the ob- 
served synchrotron emission and/or by dynamical arguments, Rawlings & Saunders 
(1991) found a nice correlation between the average power that must be supplied 
to the lobes and the power emitted by the narrow line region. Although one always 
expects some correlation between powers (they both scales with the square of the 
luminosity distance) it is the ratio of the two quantities to be interesting, being of 
order of 100. Since we also know that, on average, the total luminosity in narrow 
lines is of the order of one per cent of the ionizing luminosity, we have the remark- 
able indication that the power carried by the jet (supplying the extended regions of 
the radio-source) and the power produced by the accretion disk (illuminating the 
narrow line clouds) are of the same order. 

Celotti, Padovani and Ghisellini (1997) later confirmed this by calculating the 
kinetic power of the jet at the VLBI scale (see Celotti & Fabian 1993) and the broad 
line luminosity (assumed to reprocess ~ 10% of the ionizing luminosity). 

A possible explanation involves the magnetic field being responsible for both the 
extraction of spin energy of a rotating black hole and the extraction of gravitational 
energy of the accreting matter. Assume in fact that the main mechanism to power 
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the jet is the Blandford-Znajek (1977) process: 



Ljet * (^)V B (3i? s ) 2 c 



(1) 



where (a/m) is the specific black hole angular momentum (~ 1 for a maximally 
rotating Kerr hole), Ub is the magnetic energy density and R s is the Schwarzchild 
radius. Note that Eq. 1 has the form of a Poynting flux. Assume now that most of 
the luminosity of the accretion disk is produced at 3i? s . The corresponding radiation 
energy density is then U T = Ldi s k/(367r i? 2 c), leading to 



Therefore a magnetic field in equipartion with the radiation energy density of the 
disk would lead to Lj et ~ L^isk- 

3. MASS OUTFLOWING RATE 

We can estimate the ratio of the outflowing (in the jet) to the inflowing mass rate, 
since 



W = r/M in c 2 ; L jct - TM out c 2 ; - M out = J-^M in (3) 



If jets carry as much energy as the one produced by the accretion disk, we then 
obtain that the mass outflow rate is ~ 1% of the accreting mass rate (if r\ = 10% 
and T = 10). 

4. THE BLAZAR DIVERSITY 

BL Lac objects and Flat Spectrum Radio Quasars (FSRQ) are characterized by 
very rapid and large amplitude variability, power law spectra in restricted energy 
bands and strong 7-ray emission. These common properties justify their belonging 
to the same blazar class. However they differ in many other respects, such as the 
presence (in FSRQ) or absence (in BL Lacs) of broad emission lines, the radio to 
optical flux ratio, the relative importance of the 7-ray emission, the polarization 
degree, and the variability behavior. Within the BL Lac class, Giommi & Padovani 
(1994) have subdivided the objects according to where (i.e. at what frequency) the 
first broad (synchrotron) peak is located. Low energy peaked BL Lacs (LBL) show 
a peak in the IR-optical bands, while in High energy peaked BL Lacs (HBL) this 
is in the X-ray band (see, in this volume, the contributions of Costamante et al., 
Giommi et al., Pian et al.,Tagliaferri et al., Tavecchio & Maraschi, Wolter et al.). 

As the emission of all blazars is beamed towards us, so there must be a parent 
population of objects pointing in other directions. The parent populations of BL 
Lacs and FSRQs are believed to be FR I and more powerful FR II radio galaxies, 
respectively (see the review by Urry & Padovani 1995). The absence of broad emis- 
sion lines in BL Lacs is shared by FR I radio galaxies, whose nuclei are well visible 



W = u r {m s fc 



(2) 
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by Hubble Space Telescope observations (Chiaberge, Capetti & Celotti 1999). This 
suggests that in FR I and BL Lac objects broad emission lines are intrinsically 
weaker than in more powerful objects. 

5. THE RE - UNITED BLAZARS 

Fossati et al. (1998) found that the SED of all blazars is related to their observed lu- 
minosity. There is a rather well defined trend: low luminosity objects are HBL-like, 
and furthermore their high energy peak is in the GeV-TeV band. As the bolomet- 
ric luminosity increases, both peaks shift to lower frequencies, and the high energy 
emission is increasingly more dominating the total output. Q Ghiscllini et al. (1998), 
fitted the SED of all blazars detected in the 7-ray band for which the distance 
and some spectral information of the high energy radiation were available. They 
found a correlation between the energy 7 pe akTO e c 2 of the electrons emitting at the 
peaks of the spectrum and the amount of energy density U (both in radiation and 
in magnetic field), as measured in the comoving frame: 7 pG ak oc U~ os . This indi- 
cates that, at 7p Ca k, the radiative cooling rate 7(7 p0 ak) oc jp^U ~const. It also 
suggests that this may be due to a "universal" acceleration mechanism, which must 
be nearly independent of 7 and U: in less powerful sources with weak magnetic 
field and weak lines the radiative cooling is less severe and electrons can be accel- 
erated up to very high energies, producing a SED typical of a HBL. The paucity 
of photons produced externally to the jet leaves synchrotron self-Compton as the 
only channel to produce high energy radiation. At the other extreme, in the most 
powerful sources with strong emission lines, electrons cannot be accelerated to high 
energies because of severe cooling. Their spectrum is therefore peaked in the far IR 
and in the MeV band. In these sources the inverse Compton scattering off externally 
produced photons is the dominant cooling mechanism, producing a dominant 7-ray 
luminosity. 

5.1. Powers 

For the same sample of blazars fitted in Ghisellini et al. (1998) we can estimate 
the powers radiated and transported by jets in the form of cold protons, magnetic 
field and hot electrons and/or electron-positron pairs. Since the model allows to 
determine the bulk Lorentz factor, the dimension of the emitting region, the value 
of the magnetic field and the particle density, we can then determine 

B 2 

L p = TTR 2 T 2 /3cn'm p c 2 ; L c = ttR 2 T 2 /3cn' e (-/}m c c 2 ; L B = ttR 2 T 2 /3c — (4) 

877 

where n' p and n' e are the comoving proton and lepton densities, respectively, R is 
the cross section radius of the jet, and (j)m e c 2 is the average lepton energy. These 

X A note of caution: the limited sensitivity of EGRET (onboard CGRO) and ground based 
Cherenkov telescopes allows to detect sources which are in high states. Therefore the trend of 
more high energy dominated spectra as the total power increases strictly refers to high states. 
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powers can be compared with the radiated one estimated in the same frame (in 
which the emitting blob is seen moving). The power radiated in the entire solid angle 
is thus L T = L'j.T 2 (the same holds for the power L syn emitted by the synchrotron 
process). All these quantities are plotted in Fig. 2 (Celotti & Ghisellini 2000, in 
prep.). In this figure hatched areas correspond to BL Lac objects. Several facts are 
to be noted: 

• If the jet is made by a pure electron-positron plasma, then the associated 
kinetic power is L e . However, we note that L c <C L T posing a serious energy 
budget problem. 

• If there is a proton for each electron, the bulk kinetic power L p <~ 10L r . This 
corresponds to an efficiency of ~ 10% in converting bulk into random energy. 
The remaining 90% is therefore available to power the radio lobes, as required. 

• The power in the Poynting flux, Lb, is of the same order of L Cl indicating that 
the magnetic field is close to equipartition with the electron energy density. 
This suggests that, on these scales, the magnetic field is not a prime energy 
carrier, but is a sub-product of the process transforming bulk into random 
energy. 

6. INTERNAL SHOCKS 

The central engine may well inject energy into the jet in a discontinuous way, with 
individual shells or blobs having different masses, bulk Lorentz factors and energies. 
If this occurs there will be collisions between shells, with a faster shell catching up 
a slower one. This idea has become the leading model to explain the emission of 
gamma-ray bursts, but it was born in the AGN field, due to Rees (1978) (see also 
Sikora 1994). 

• Location — The 7-ray emission of blazars and its rapid variability imply 
that there must be a preferred location where dissipation of the bulk motion 
energy occurs. If it were at the base of the jet, and hence close to the accretion 
disk, the produced 7-rays would be inevitably absorbed by photon-photon 
collisions, with associated copious pair production, reprocessing the original 
power from the 7-ray to the X-ray part of the spectrum (contrary to obser- 
vations). If it were far away, in a large region of the jet, it becomes difficult to 
explain the observed fast variability, even accounting for the time-shortening 
due to the Doppler effect. The region where the radiation is produced is then 
most likely located at a few hundreds of Schwarzchild radii (~ 10 17 cm) from 
the base of the jet, within the broad line region (see Ghisellini & Madau 
1996 for more details). The extra seed photons provided by emission lines 
enhance the efficiency of the Compton process responsible for the 7-ray emis- 
sion. This is indeed the typical distance at which two shells, initially separated 
by i?o ~ 10 15 cm (comparable to a few Schwarzschild radii) and moving with 
r - 10 and T ~ 20 would collide. 
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FIGURE 1. Histograms of the powers carried by the jet in protons, total radiation, 
synchrotron radiation, magnetic field and relativistic electrons, from top to bottom. 
Hatched areas correspond to BL Lac objects. The electron distribution was assumed 
to extend down to 7 m i n ~ 1. From Celotti & Ghisellini (2000, in prep.) 
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FIGURE 2. Cartoon illustrating the internal shock scenario. The intermittent ac- 
tivity of the central engine produces two shells, initially separated by i?o- The faster 
one will catch up the slower one at R ~ T 2 i?o- 



• Variability timescales — In fact if the initial separation of the two shells 
is i?o and if they have Lorentz factors Ti, T2, they will collide at 

or 2 

If the shell widths are of the same order of their initial separation the time 
needed to cross each other is of the order of R/c. The observer at a viewing 
angle 9 ~ 1/T will see this time Doppler contracted by the factor (1— (3 cos 9) ~ 
r~ 2 . The typical variability timescale is therefore of the same order of the 
initial shell separation. If the mechanism powering GRB and blazar emission 
is the same, we should expect a similar light curve from both systems, but 
with times appropriately scaled by the different Ro, i.e. the different masses 
of the involved black holes. 

• Efficiencies — As most of the power transported by the jet must reach the 
radio lobes, only a small fraction can be radiatively dissipated. The efficiency 
77 of two blobs/shells for converting ordered into random energy depends on 
their masses mi, and bulk Lorentz factors T\, T%, as 

mi+m 2 

7? = 1 — I t (6) 

' ; Timi + r 2 TO 2 
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where I 1 / = (1 — f3j) 1 I 2 is the bulk Lorentz factor after the interaction and 
is given by (see e.g. Lazzati, Ghiscllini & Celotti 1999) 



The above relations imply, for shells of equal masses and T 2 = 2Li = 20, 
T f = 14.15 and r) = 5.7%. 

Efficiencies X] around 5-10% are just what needed for blazar jets. 

• Peak energies? — In the rest frame of the fast shell, the bulk kinetic en- 
ergy of each proton of the slower shell is ~ (V — l)m p c 2 , where L' ~ 2. This 
is what can be transformed into random energy. Assume now that the elec- 
trons share this available energy (through an unspecified acceleration mecha- 
nism). In the comoving frame, the acceleration rate can be written as Eheat ~ 
(r' — l)mj,c 2 jt'heaf ^ ne typical heating timescale may correspond to the time 
needed for the two shells to cross, i.e. t' heat ~ AR' /c ~ R/(cT), where AR' is 
the shell width (measured in the same frame). The heating and the radiative 
cooling rates will balance for some value of the random electron Lorentz factor 



The agreement of the above simple relation with what can be derived from 
model fitting the SED of blazars is surprisingly good (see Ghiscllini 2000). 

• Radio flares — Collisions between shells may (and should) happen in a 
hierarchical way. As an illustrative example, assume that one pair of shells 
after the collision moves with a final Lorentz factor Ti = 14 (this number 
corresponds to T = 10 and 20 for the two shells before the interaction). The 
collision produces a flare -say- in the optical and 7-ray bands. After some 
observed time At two other shells collide and another flare is produced. As- 
sume that the final Lorentz factor is now r 2 = 17 (corresponding to an initial 
L = 10 and 30 before collision). Since the second pair is faster, it will catch 
up the first one after a distance (from eq. 5) R ~ 1200cAt. A time separation 
of At ~ a day between the two flares then corresponds to R ~ 1 pc, i.e. the 
region of the radio emission of the core. Due again to Doppler contraction, 
this radio flare will be observed ony a few days after the second optical flare. 
Since the ratio T2/T1 is small, the efficiency is also small (at least a factor 
10 smaller than the firsts shocks). There is then the intriguing possibility of 
explaining the birth of radio blobs after intense activity (i.e. more than one 
flare) of the higher energy flux. Radio light-curves should have some memory 
of what has happened days-weeks earlier at higher frequencies. 



gilVni + /3 2 r 2 m 2 
T\m\ + T 2 m 2 



(7) 
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7. CONCLUSIONS 

Here I will dare to assemble different pieces of information gathered in recent years 
in a coherent, albeit still preliminary, picture. 

There is a link between the extraction of gravitational energy in an accretion disk 
and the formation and acceleration of jets, since both have the same power. Objects 
of low luminosity accretion disks also lack strong emission lines, suggesting that it 
is the paucity of ionizing photons, not of gas, the reason for the lack of strong lines 
in BL Lacs. Correspondingly, this implies that, if FR I are the parents of BL Lacs, 
they also have intrinsically weak line emission (i.e. no need for an obscuring torus). 
Despite the fact that the jet power in blazars spans at least four orders of magnitude, 
the average bulk Lorentz factor is almost the same, suggesting a link between the 
power and the mass outflowing rate: their ratio is constant. In the region where 
most of the radiation is produced, the jet is heavy, in the sense that protons carry 
most of the bulk kinetic energy. There the jet dissipates ~ 10% of its power and 
produces beamed radiation. The power dissipated at larger distances is much less, 
and therefore the jet can transport ~ 90% of its original power to the radio extended 
regions. One way to achieve this is through internal shocks, which can explain why 
the major dissipation occurs at a few hundreds Schwarzchild radii, why the efficiency 
is of the order of 10%, and give clues on the observed variability timescales and 
even on why electrons are accelerated at a preferred energy. The spectral energy 
distribution of blazars depends on where shell-shell collisions take place, and on 
the amount of seed photons present there. Even in a single source it is possible that 
the separation of two consecutive shells is sometimes large, resulting in a collision 
occuring outside the broad line region. In this case the corresponding spectrum 
should be produced by the synchrotron self-Compton process only, without the 
contribution of external photons: we then expect a simultaneous optical-7-ray flare 
of roughly equal powers (but with the self-Compton flux varying quadratically, see 
Ghiscllini & Maraschi 1966). This is what should always happen in lincless BL Lac 
objects. On the other hand, if the initial separation of the two shells is small (or 
the T-factor of the slower one is small), the collision takes place close to the disk. 
X-rays produced by the disk would then absorb all the produced 7-rays and a pair 
cascade would develop, reprocessing the power originally in the 7-ray band mainly 
into the X-ray band. We should therefore see an X-ray flare without accompanying 
emission above Fm c c 2 . 

Pairs of shells which have already collided can interact again between themselves, 
at distances appropriate for the radio emission. This offers the interesting possibility 
to explain why the radio luminosity is related with the 7-ray one, and why radio 
flares are associated with flares at higher frequencies. Work is in progress in order 
to quantitatively test this idea against observations. 
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